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ABSTRACT 

In the lenticular galaxy NGC 1023 a new population of star clusters ("faint 
fuzzies") was recently discovered by Larsen & Brodie. These clusters are found 
inside the disc and are faint (23 < V < 24 mag) and extended with effective radii 
of r cff 7 to 15 pc. 

We present here N-body calculations of a likely formation-scenario through merg- 
ing star clusters in clusters of star clusters (super-clusters). Such super-clusters 
are observed to form in interacting galaxies. 

The resulting merger objects have masses comparable to the "faint fuzzies" and 
show large effective radii (r eff > 7 pc). Even though these objects are suffering 
from strong tidal forces they are able to survive and reach the estimated ages of 
the extended star clusters in NGC 1023. 

Subject headings: galaxies: star clusters - galaxies: interaction - galaxies: 
individuals (NGC 1023) - methods: N-body simulations 

1. Introduction 

During a recent search for globular clusters (GC) in NGC 1023, a lenticular 
galaxy at approximately 9 Mpc distance, Larsen & Brodie (2000) found faint objects 
(23 < V < 24 mag) that have GC luminosities but effective radii (= half-light radii) of 7 to 
15 pc which is much larger than ordinary GCs that have effective radii of 2 to 3 pc. Using 
spectra obtained with the Keck telescope they (Brodie & Larsen 2002) confirmed that these 
objects belong to NGC 1023 and have ages of about > 7 to 8 Gyr, are moderately metal rich 
([Fe/H] = —0.58 ± 0.24) while showing a super-solar a-element abundance ([a/Fe] = +0.3 
to +0.6). While GC are usually connected with the spherical component of a galaxy, these 
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new objects appear to be situated in the disc of NGC 1023. The distribution of these 
objects shows rotation around the galaxy. The authors conclude that they have found a 
new class of star clusters which they call "faint fuzzies" . 

NGC 1023 has a dwarf companion (NGC 1023A) with which it is steadily interacting, 
and due to its lenticular appearance it probably suffered from interactions and mergers in 
the past leading to strong star formation events. 

Interactions between gas-rich galaxies lead to strong star formation events. For exam- 
ple, HST-images of the Antennae reveal (Whitmore et al. 1999) that the knots of intense 
star formation produce clusters of massive young star clusters. These "super-clusters" (= 
clusters of star clusters; not to be confused with super stellar clusters (SSC), which are 
individual massive star clusters) appear to contain dozens to hundreds of massive star 
clusters within a region spanning a few hundred pc to a kpc in diameter. 

Here we report N-body results concerning the dynamical evolution of such super- 
clusters. The resulting merger objects have similar properties as the "faint fuzzies". 

2. Simulations 

We use the particle-mesh code Superbox (Fellhauer et al. 2000) which incorporates a 
hierarchical grid architecture allowing high resolution at the places of interest. 

We model the single star clusters as Plummer spheres (Plummer 1911, numerical 
realisation: Aarseth et al. 1974) with a Plummer radius of 4 pc, which corresponds to 
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the mean half-light radius found for the individual young star clusters in the Antennae 
(Whitmore et al. 1999). iVo = 262 star clusters with a mass range of 10 4 — 10 6 M Q following 
a power-law mass spectrum, 



n(M cl ) oc M^ 1 , (1) 

are placed in a Plummer distribution with Plummer radius r^j representing the super-cluster. 
The super-cluster orbits on an eccentric orbit around the host galaxy represented by an 
analytical potential, which consists of a disc modelled as a Plummer-Kuzmin potential and 
a spherical halo component modelled as a logarithmic potential: 



$gal = $disc + $halo (2) 



GMdisc 1 2 i / D 2 



^R 2 + {a + Vz 2 + b 2 ) 2 



with M disc = 10 11 M , a = 3 kpc, b = 0.3 kpc, v = 200 km/s and i? ga i = 50 kpc which 
sums up to an almost flat rotation curve with a rotation speed of 220 kms" 1 . 



Several simulations with different values for the scale length r?J of the super-cluster, dif- 
ferent initial masses M sc of the super-cluster and different orbits around the host galaxy are 
performed to study the evolution of the resulting merger object. The computational setup 
closely follows that used in previous studies (Fellhauer et al. 2002; Fellhauer & Kroupa 2002). 



It is possible to follow the evolution with a particle-mesh code that neglects dynamical 
effects of two-body relaxation, because the half-mass (bulk) two-body relaxation times of 
the single star clusters, which can be estimated by (Binney & Tremaine 1987) 

664 / M cl \ 1/2 /1M \ / ro^ 3/2 



w = hK5n){wm^) l^rjfe} Myr ' (3) 



- 5 - 



is 800 Myr for a 10 4 M star cluster ranging up to 4.4 Gyr for a 10 6 M star cluster, 
while the merging timescale is much shorter (Fellhauer et al. 2002). Furthermore, as shown 
below, the resulting merger objects have relaxation times of a Hubble-time or longer. 

3. Results 

In a previous paper (Fellhauer et al. 2002) it was shown that the star clusters inside 
a super-cluster merge within a few crossing times of the super-cluster (typically 50 to 
500 Myr). The resulting merger objects are stable and compact. Their masses range 
from being similar to GCs, in cases where the tidal field removes a large fraction of 
the individual clusters leaving merger objects consisting of only a few star clusters. On 
the other hand, when the super-cluster is sufficiently concentrated and the tidal field 
is correspondingly weak the merger object evolves to become an ultra-compact dwarf 
galaxy (UCD) (Fellhauer & Kroupa 2002), a new class of objects discovered recently in 
the Fornax galaxy cluster (Phillipps et al. 2002; Hilker et al. 1999). To show the different 
results of a wide range of our simulations we plot the effective radii and the masses of 
the merger objects as a function of the Plummer radius of the super-cluster (Fig. 1). 
Even if tidal forces are strong, compact and massive super-clusters evolve into extended 
objects with lifetimes significantly longer than 5 Gyr although they orbit well inside the disc. 

Here we give the results of a simulation (OC03, filled box in Fig. 1) where the 
super-cluster was placed on an orbit within the galactic disc which accounts for the scenario 
in NGC 1023. The resulting merger object is far too massive to be one of the "faint fuzzies" 
found in NGC 1023, rather resembling an object like u-Cen in our Milky Way but with two 
additional simulations we show that these results still hold for less massive objects like the 
faint extended star clusters in NGC 1023. 
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The super-cluster has initially a Plummer radius of rj^ = 50 pc and a cut-off radius of 
250 pc. It contains 262 star clusters following the power law of Eq. 1 with a total mass of 
10 7 M . The crossing time of the super-cluster is t^. = 10.4 Myr and the internal velocity 
dispersion of the clusters in the super-cluster is <r sc = 15.9 kms -1 . Additionally we choose 
the sense of rotation of all clusters in the super-cluster to be the same to investigate the 
resulting rotation-law of the merger object. A super-cluster is expected to rotate if it forms 
from a contracting and locally differentially rotating inner tidal arm. The super-cluster is 
placed on an eccentric orbit with perigalacticon at 2.1 kpc and apogalacticon at 7.5 kpc. 
The orbit is inclined such that the maximum Z-distance from the disc plane is about 2 kpc. 
The parameters are chosen to be representative of the knots seen to contain many star 
clusters in the Antennae galaxies (Fellhauer & Kroupa 2002), while the orbital inclination 
is motivated by the orbit of cj-Cen (Dinescu et al. 1999). A possible link between u-Cen 
and merged super-clusters has already been pointed out by Fellhauer & Kroupa (2002). 

Within the first 250 Myr a merger object forms with a bound mass of about 3 • 10 6 M Q . 
The half-mass radius of this object is 18 pc and the 90% Lagrangian-radius is of the order 
of 110 pc at that time. The surface density profile of the merger object follows a shallow 
power law of £ oc r~ 2 5 as shown in Fig. 2a. The 3-dimensional central velocity dispersion 
is °"3£> = 21.8 kms~ x (Fig. 3a). It also shows a maximum rotation velocity of about 5 kms -1 
at a distance of about 12 pc from the centre (Fig. 3a). The bulk relaxation time (Eq. 3) of 
the merger object is ~ 50 Gyr which is longer than the age of the universe. 

The tidal radius calculated using the equation 7-84 from Binney & Tremaine (1987) 




(4) 
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where D denotes the distance to centre of the galaxy and M gal (D) is the enclosed mass of 
the galaxy at that distance and M mo is the bound mass of the merger object, is about 84 pc 
at perigalacticon. Therefore the object becomes tidally shaped and looses mass with every 
perigalacticon passage. After 5 Gyr of simulation time the merger object has lost about half 
its mass (Fig. 5a). The central surface density drops down to 7500 M pc~ 2 . The profile can 
be fitted with a King profile with King radius (= effective radius) of 6 pc. An alternative 
fit is an exponential profile in the inner part (£ oc e _r//6pc ) and a steep power law with 
E oc r~ 5 5 for the outer part (Fig. 2b). The half-mass radius has shrunk down to 11.6 pc, 
the 90% Lagrangian-radius to 27 pc. The tidal radius is determined to be 68 pc at the last 
perigalacticon and 127 pc at the actual position. Also the central velocity dispersion and 
the maximum rotation speed has dropped as shown in Fig. 3b. The bulk relaxation time 
(following Eq. 3) of the merger object is now t rclax ps 20 Gyr. 

The two additional simulations have N = 20 star clusters initially with masses of 
10 4 M (M sc = 2 x 10 5 M ) (OC07) and 5 x 10 4 M (M sc = 10 6 M ) (OC08) respectively. 
The Plummer-radii of the super-clusters are chosen to be 20 pc. The orbits of the 
super-clusters have the same eccentricity as in the first simulation but are placed in the 
X-F-plane of the galaxy instead of being inclined. 

In the OC07 simulation a merger object containing 14 out of the initial 20 star clusters 
develops after 100 Myr. After 5 Gyr this merger object has a bound mass of 1.8 • 10 4 M Q 
(Fig. 5a) and a tidal radius of 24 pc. It's surface density profile can be fitted with a King 
profile with a central surface density of 240 M pc~ 2 and a King radius (which corresponds 
to the effective radius) of 5.9 pc. Fitting an exponential profile to the inner part gives 
E exp = 250 M pc -2 and an exponential scale length of about 4 pc. The outer part can 
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again be fitted with a steep power law Eocr 5 - 9 (Fig. 4a). The central line-of-sight velocity 
dispersion is 1.4 kms -1 . 

The OC08 simulation shows a merger object which developed out of 18 star clusters. 
Again the merging process is finished after about 100 Myr. The mass of the merger object 
after 5 Gyr is 4 • 10 5 M (Fig. 5a) and the tidal radius is 70 pc. Fitting a King profile 
to the surface density distribution (Fig. 4b) gives the following values: central surface 
density of 1500 M pc -2 and an effective radius of 7.5 pc. The exponential fit has the values 
1400 M Q pc~ 2 and r cxp = 7.5 pc. The power law fit to the outer part has an power-law index 
of -6.0. 

The above results (total mass and central surface density) can be converted to total 
blue luminosity and central surface brightness using mass-to-light ratios taken from a single 
stellar population model of Starburst99 (Leitherer et al. 1999). The programme adopts 
a Salpeter IMF. Therefore, we choose a minimum stellar mass of 0.3 M and a maximum 
mass of 100 M . Note that we truncate the Salpeter IMF at 0.3 M Q . This roughly accounts 
for the flattening of the IMF below 0.5 M (Kroupa 2002). The evolution of M/Lb with 
time is shown in Fig. 5b as obtained as a result of stellar evolution from Starburst99. 
We stop the photometric evolution from Starburst99 when M/L B = 3.0 is reached 
because Starburst99 overestimates the M/L at high ages quite substantially. Instead we 
take a linear fading with time as proposed by the models of Bruzual & Chariot (1993) of 2 
magnitudes from 3 to 15 Gyr. This enables us to place the evolutionary track of our merger 
objects into a central surface brightness, //o,b, vs absolute photometric B-band magnitude, 
Mb, diagram (the Kormendy diagram; Fig. 6). The endpoint of the evolution of our merger 
objects is also the region where the "faint fuzzies" are located. The faint fuzzies have total 
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luminosities of about M v = —7 (Larsen & Brodie 2002). We estimate their location in 
the Kormendy diagram in the following way. First we set My + 1.0 = M# (i.e. adopting 
(B — V) — 1.0 as an admissible approximation for the diagram) then we convert the total 
luminosity into a total mass using M/L = 3.0. With this mass we calculate the tidal radius 
in our model potential and orbit (r t id a i ~ 50 pc) and adopt an effective (or King-) radius of 
10 pc. Using a King profile we estimate the central surface brightness to be // ,b ~ 22.3 mag. 

4. Conclusions 

SO (lenticular) galaxies are commonly believed to have suffered from interactions 
with other galaxies which transformed their morphology from a late-type spiral into their 
early-type form (Kennicutt 1998; Abraham & van den Bergh 2001). Therefore the type 
of NGC 1023 tells us that this galaxy suffered maybe several but at least one strong 
interaction (merger) with another galaxy. As is seen in interacting galaxies these events lead 
to an enhanced star formation activity and to an enhanced rate of super nova explosions 
of massive stars which accounts for the super-solar enhancement of a-elements. Many of 
the new stars mostly form in star clusters which are also grouped in super-clusters, as 
is demonstrated by the Antennae and Stephan's Quintet (Gallagher et al. 2001). These 
super-clusters are bound entities and form, through successive mergers of their constituent 
star clusters, new objects which are larger than ordinary GCs. Our simulations show 
that even strong tidal fields are not able to destroy these objects and they survive for a 
Hubble-time. 

The "faint fuzzies" found in NGC 1023 may have therefore formed as super-clusters 
during the interactions and the resulting star bursts that shaped NGC 1023. Like in 
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the Antennae the super-clusters were born within the galaxy and not in the outer tidal 
tails. Therefore the "faint fuzzies" , which one can call extended globular clusters or small 
ultra-compact dwarf galaxies by virtue of their long two-body relaxation time, stay within 
the optically visible galaxy and show a similar rotation pattern like the disc. These objects 
can be old because they are able to survive until the present, they may be rotating and 
have different populations of stars if old field stars were captured by the young super-cluster. 

Given the ubiquity of young super-clusters in the Antennae and Stephan's Quintet, 
it appears reasonable to expect that such objects form whenever strong star bursts occur. 
It is therefore exciting, but retrospectively not very surprising, that "faint fuzzies" were 
discovered in a lenticular galaxy. 

Super-clusters that form in outer tidal tails may leave the mother galaxy either by 
being ejected by the merging galaxies or become unbound due to a tidal field from the 
galaxy cluster. Such merger objects will appear similar to the UCDs in Fornax (Fig. 7). A 
related object may be u-Cen in our Milky Way, which is quite different and more massive 
than the other GCs. It incorporates at least two different populations of stars and shows 
a rotation-law resembling that of our merger object (Freeman 2001). oCen is quite old 
and could therefore be a relic of a star-burst during the formation of the Milky Way. We 
are investigating the hypothesis that u-Cen may also be related to super-clusters. Fig. 7 
summarises in a schematic way the possible evolutionary paths and resulting objects that 
may form under this hypothesis. 

MF acknowledges financial support through DFG-grant FE564/1-1. MF also 
acknowledges the support of the European Commission through grant number HPRI-1999- 
CT-00026 (the TRACS programme at EPCC). 
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Fig. 1. — Effective radii, r e a, (a) and masses, M mo , (b) of our merger objects evaluated at 
5 Gyr vs. Plummer radius of the initial super-cluster r^j. Crosses: large initial super-cluster 
mass (M sc > 10 7 M Q ) and weak tidal field (distance at perigalacticon _Rp Cri > 10 kpc); boxes: 
large mass and strong tidal field (distance at apogalacticon R apo < 10 kpc); triangles: low 
mass (M sc < 10 7 M ) and strong tidal field; tri-pointed stars: low mass and weak tidal field. 
Filled symbols show the models discussed in the paper. 
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Fig. 2. — Projected surface density profile of the merger object. The vertical lines are as 
follows: dashed: half-mass and 90%-Lagrangian radius; solid: tidal radius at the position of 
the merger object (r t ) and at the last perigalacticon (r t:P ). (a): Measured at t — 350 Myr 
(shortly after formation of the merger object). The dashed line is a power law with E oc r -2 5 . 
The adopted mass-to-light ratio to calculate the surface-brightness is 0.3. (b): Profile at 
t — 5 Gyr. Dot-dashed line shows the fitted King profile, solid line is the exponential fit 
for the inner part and dashed line is the fitted power law for the outer part. The adopted 
mass-to-light ratio is 3. 
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Fig. 3. — Velocity distribution in the merger object; circles show the 3-dimensional velocity 
dispersion computed in concentric shells around the centre; three-, four- and six-pointed 
stars are the projected line-of-sight velocity dispersions along the x-, y- and z-coordinate 
axis evaluated in concentric rings. Squares show the mean rotational velocity again measured 
in concentric shells. Note the similarities to the rotation-law observed for uj-Cen (Freeman 
2001). Vertical line shows the tidal radius at perigalacticon. (a): Measured at t — 350 Myr 
(shortly after formation of the merger object), (b): Distributions at t — 5 Gyr. 
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Fig. 4. — Projected surface density profile of the merger objects. The vertical lines are the 
calculated tidal radii at perigalacticon (r tiP ) and at the actual distance (r t ) (at t — 5 Gyr) of 
the merger objects. The solid curves show the fitted exponential profiles, the short dashed 
the fitted King profile and the long dashed the fitted power law with the values according 
to the main text, (a) OC07: Simulation with M sc = 2 • 10 5 M ; (b) OC08: Simulation with 
M sc = 10 6 M . 
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Fig. 5. — (a): Mass of the merger objects vs. time. From top to bottom: M sc = 10 7 M 
(OC03), M sc = 10 6 M (OC08) and M sc = 2 • 10 5 M (OC07). (b): Evolution of the 
mass-to-light ratio M/Lb taken from Starburst99 (dashed line) and the "frozen" value 
(M/Lb = 3.0) used to calculate the luminosities (solid line). 



-18- 




Fig. 6. — Central surface brightness, //o,b, vs absolute photometric B-band magnitude, M B 
(the Kormendy diagram). The stars are data for Milky- Way globular clusters (GC, Har- 
ris 1996), Local-Group dwarf galaxies (Mateo 1998) are displayed as filled boxes, elliptical 
galaxies (E, Peletier et al. 1990) are shown as open boxes. The positions of disk galaxies 
are indicated by "Spirals" (cf. Ferguson & Binggeli 1994). The newly discovered UCDs 
are shown by arrows with lower-limits on //b (Phillipps et al. 2001). Three "knots", (KnS, 
Kn430, Kn225, tables 1 and 2 in Whitmore et al. 1999), observed in the interacting An- 
tennae galaxies are shown as stars. Knots KnS and Kn430 are roughly 10 Myr old so that 
M/L B w 0.035, while Kn225 is about 500 Myr old (M/L B « 0.476). The open six-pointed 
star symbol shows the merged super-clusters predicted by Kroupa (1998) to result from the 
dynamical evolution of such knots, which consist of dozens to hundreds of young star clus- 
ters. The large fuzzy symbol shows the approximate location of the faint fuzzy star clusters. 
The lines show the evolution of our merger objects from the time of formation until 5 Gyr 
using M/L B from Fig. 5b and an additional linear fading for ages > 3 Gyr (from top-left to 
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Fig. 7. — Schematic of possible evolutionary tracks of super-clusters. Question marks signify 
on-going research. 



